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A direct connection between the magnitude of the violation of the fluctuation-response relation 
(FRR) and the rate of energy dissipation is presented in terms of field variables of nonequilibrium 
systems. Here, we consider the density field of a colloidal suspension either in a relaxation process 
or in a nonequilibrium steady state driven by an external field. Using a path-integral representation 
of the temporal evolution of the density field, we find an equality that relates the magnitude of 
the violation of the FRR for scalar and vector potentials of the velocity field to the rate of energy 
dissipation for the entire system. Our result demonstrates that the violation of the FRR for field 
variables captures the entropic component of the dissipated free energy. 
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In recent years, considerable effort has been made to 
understand the nature of fluctuations in nonequilibrium 
systems. Of particular note are studies on the fluctu- 
ation theorems [II, I3j S] j the Jarzynski equality ^] , the 
Hatano-Sasa equality [5| , and the other related equalities 
that hold exactly for systems driven far from equilibrium. 
In addition to these works on rather established topics, 
further novel results regarding the properties of nonequi- 
librium systems are still being obtained 7]. 

With regard to fluctuations in nonequilibrium systems, 
we recently reported on the physical significance of the 
violation of the fluctuation-response relation (FRR) in 
systems far from equilibrium Q . The FRR is a rela- 
tion between the response function and the correlation 
function of gross variables. This relation has been firmly 
established for systems near equilibrium [l^], and it is 
known to be violated in various nonequilibrium systems 
[1, B ■ In previous works d, 0] , we presented sev- 

eral equalities that directly relate the amount by which 
the FRR is violated and the rate of energy dissipation for 
systems described by Langevin equations. These equali- 
ties have been extended in several ways 13l . |14| and have 
been experimentally verified 15, 3, 17 1. In these studies, 
violation of the FRR is considered in terms of mechani- 
cal degrees of freedom. It has not yet been demonstrated 
whether a similar equality also holds for field variables, 
e.g., the density field of the particles. In the present 
Rapid Communication, we demonstrate that in fact vio- 
lation of the FRR for field variables is also directly con- 
nected to the dissipation rate. Our result provides new 
insight into violation of the FRR in nonequilibrium sys- 
tems as well as the proposal of a novel method to experi- 
mentally probe energy dissipation in nonequilibrium soft 
materials. 

For the sake of concreteness, here we consider n spher- 
ical colloidal particles suspended in an aqueous solution 
of temperature T and confined to a three-dimensional 



box [18|. The particles are subject to a potential field, 
U{x), and a static driving force, F{x). We suppose that 
F{x) cannot be written in a gradient form, and it thus 
represents a nonconservative force. The particles are in- 
teracting with each other via a potential represented by 
W{x). For simplicity, we assume W{x) — W{—x). For 
such a system, the motion of the i-th particle is described 
by the following Langevin equation: 

n 

jx, = Fix,) - S/U{x,) - J2 ^Wix, - Xj) + I,. (1) 

Here, 7 denotes the friction coefficient of a particle, and 
^j(t) represents zero- mean white Gaussian noise satisfy- 
ing {^i{t)$j{s)') = 2'jTSij\6{t — s), where I is unit matrix. 
The Boltzmann constant is set to unity. For simplicity, 
we do not consider the hydrodynamic interactions among 
particles in the present Rapid Communication. Note that 
the model ([1]) has been considered in many contexts, in- 
cluding the rheology of colloidal suspensions [l^, glass 
transitions |2Q] . and fundamental studies of nonequilib- 



rium steady states [21 



Here, we are interested in fluctuations and the response 
of the density held of the particles denoted by p{x, t). In 
order to study the response of the system, one can in- 
troduce a perturbative external force, which we write as 
f{xi,t) = -e\/(j){xi,t) + AV X A{xi,t) (e, A < 1), and 
add this to the r.h.s. of In particular, two types of 
responses are important in the following. One is the re- 
sponse of a scalar field, X[[/9]](cc, t), which is a functional 
of p, to the scalar field e(j), written in the form 

{Xix,t))^, = {X{x,t)), 

-e I dt' I dxRx{x,t-x',t')<j){x',t') + 0{e^),{2) 

which defines a response function Rx- Here, (■)^ ^ ((■)o) 
represents the ensemble average in the presence (absence) 



2 



of the perturbation fields, and denotes a sufficiently 
large region that contains the box. We suppose that 
and A vanish at the boundary of O. The other type of 
response is that of a vector field, Y'[[p]](a;, t), which is 
also a functional of p, to the vector potential \A. This 
response is written in the form 

(r(a;,<))o_, = (y(a;,i))o 

-A I dt' [ dxQY{x,t;x',t')-A{x',t') + 0{\^),iS) 
Jo Jn 

which defines a response matrix Qy. Note that 
Rx{x,t;x' ,t') and QY{x,t;x' ,t') are defined only for 
t > t', while we set them to zero for t < t' , in accor- 
dance with causality. 

Next, fluctuations of the density field can be charac- 
terized by the correlation function. Representing scalar 
functionals oi phy X[[p]]{x,t) and X'[[p]]{x' ,t'), the cor- 
relation function is defined as 

Cx,x'{x,t;x',t') = {X{x,t)X'{x',t'))o- 

Then, for vector functionals of p, denoted by l^[[/3]](a;, t) 
and Y'[[p]]{x',t'), we define 

DY^Y'{x,t-x\t') = {Y{x,t)Y'ix',t'))^ 

as the correlation matrix. 

In the present Rapid Communication, we show that 
the quantities defined above are directly connected to 
the total rate of free energy dissipation in the system, 
defined as 



dt ./o 



dxj{x,t)-F{x), (4) 



where H[p{-,t)] represents the free energy of the system 
defined in and j{x, t) is the mass flux defined in PH)) . 
Our central result is the following [2^ : 



{Jn{t))Q^j / dxlCv-v j{x,t;x,t) - TRv{x,t;x,t) 
+tr [Dw.-vxj (x, t; x, t) - TQw{x, t; x, t)] }. (5) 

Here, V and W represent scalar and vector potentials of 
the velocity field, j{x,t)/p{x,t), defined as 

dx'g{x;x')V-(^^^), (6) 



(7) 



V[[p]]{x,t) ^ 

W[[p]]ix,t) = - [ dx'gix;x')V' X 
Jn 



p{x',t)j 
3{x',t) 
p{x',t) 



where g[x] x') is the Green function of the Laplacian op- 
erator in O, defined as ^g{x\ x') — —S{x — x'), with the 
boundary condition g{x; x') = for a; G 9^1. 

Interestingly, the r.h.s. of (O can be interpreted as the 
magnitude of the violation of the following FRRs: 



Cv-v.j{x,t;x',t') = TRv{x,t;x',t' 



(8) 



As demonstrated below, ([8]) and ([9]) hold when the sys- 
tem is in equilibrium. However, it has been argued that 
when the system is far from equilibrium, the FRRs ^ 
and (O do not hold in general d, [^, [ill, \lM ■ The equality 
([5]) demonstrates that there is a direct relation between 
the rate of energy dissipation and the magnitude of the 
violation of the FRR for V and W for a system far from 
equilibrium. 

It is also noteworthy that the r.h.s. of the equality ^ 
consists entirely of experimentally measurable quantities. 
Thus, using this equation, we can investigate the ener- 
getic aspects of a colloidal suspension driven far from 
equilibrium without the need for a detailed study of the 
force fields such as F,U, and W. In fact, all the quan- 
tities on the r.h.s. of (O, except for 7, can be deter- 
mined through measurements of the density field and 
the flux by considering cases with and without the per- 
turbation fields. For a colloidal suspension, it might be 
possible to make such measurements by using a confocal 
microscopy [111 , for example, and precisely calibrated ex- 
ternal fields can be created using a variety of techniques, 
such as those employing optical trap arrays [2^. Fur- 
thermore, the value of 7 can be estimated with various 
well-established methods, e.g., using the Einstein-Stokes 
relation or through measurement of the high-frequency 
response [Ij. 

We now derive ([5]). We start with the fluctuating hy- 
drodynamic equation describing the time evolution of the 
density field, pix,t), which has been obtained by several 



authors [13, IM [3, H 



d_ 

dt 



p{x,t) 
j{x,t) 



-^■3{x,t), 



(10) 



p{x,t) 



Fix)-V^-^^l^-eVc^{x,t) 



+XV X A{x,t) 



5p{x,t) 
\/p{x,t) 
1 



oC{x,t),{ll) 



where 7i represents the free energy functional, defined as 

n[if\ = I dxU{x)ip{x) +T dxip{x)[liiip{x) ~1] 
Jn Jn 

+ 1 [ dx [ dy^{x)W{x-y)^{y), (12) 
^ Jn Jn 

and C,{x^t) represents white Gaussian noise satisfying 
{C[x, t)) = and {C{x, t)C(y, s)) = 2^T\Six - y)5{t - s). 
The multiplication o is interpreted according to the Ito 
rule for both space and time [28]. Next, using the Green 
function, g{x;x'), and the Helmholtz theorem, (fTOll and 
PT|) can be rearranged as 



, , ^5H[p(-,t)] ^, , 
ri{x,t) = V ' ~ F{x) 



7 



5p{x,t) 



p{x,t) 



V X (t{x, t) 



7 



dx'\7g{x;x')-^p{x',t) 



^v^j{x,t;x',t') = TQwix,t;x',t'), (9) 



Pix,t) JO 

eV(t){x,t) - AV X A{x,t), 



(13) 



3 



where q{x, t) = p{x, o ^(a;, t), and (t{x, t) is a vec- 

tor field satisfying V ■ <T{x,t) = 0. Thus, (T{x,t) can 
be specified by two independent fields, which we write 
ai{x,t) and a2{x,t). Let p denote {p,<7i,a2) collec- 
tively. It can be shown that the noise t]{x, t) is zero- 
mean white Gaussian noise satisfying {r]{x,t)T]{x' , t')) — 
2'yT\5{x - x')6(t - t')/p{x, t). The relation ^ defines 
a transformation between ri{x,t) and p{x,t). Note that 
the functional derivative Sr)/6p is independent of e(j) and 
AA, because these fields are independent of p. Therefore, 
the Jacobian of the transformation, J[[p]] = det{Sri/Sp), 
is also independent of the perturbative fields. We thus 
obtain the path-integral representation for the transition 
probability of p from t = to t = t under a given initial 
condition, pg = p(-,0), as 

PeMiipMPa]) - J[[pm[p]]eM-PS,A[p]]), (14) 

where N[[p]] = llxIlt[iPi^^i)/^'^^'Pfd.^dt]^^'^ is a nor- 
malization constant for the path probability of t], and 
P = T~^. Note that iV[[p]] is also independent of the 
perturbation fields. The action, S'e_A[[p]], is given by 



Se.xM] = ^ I dt ( dxp{x,t) 



v«:41-F(.) 



pix,t) 



Sp{x,t) 

V X cr{x,t) + / dx'Vg{x;x')—p{x\t) 



+£V(f){x,t) - AV X A{x,t) 



Our next task is the calculation of the average of the 
functionals X[[/9]](a;, <) and l^[[p]](a;, t). The average of 
X is defined as 



VpX[[p]]ix,t)P,[p,]P,,xi[[pmPo])- 



The initial distribution for the density fields, Pi[po]j can 
be chosen arbitrarily at this point. We then study sep- 
arately the effects of the two types of perturbations, e(j) 
and \A. First, we consider the former by setting A = 0. 
After setting t = r/2, a straightforward calculation leads 
to 

{X{x,r/2))^„^{X{x,r/2)), 

-e f dt' f dx'^[Cx.v-pG{x,T/2;x',t') 
Jo Jn ^7 

+jCx.-v-j{x, t/2; x', t')] 0(a;', t') + 0{e^), (15) 

where Gix,t) = F{x) - V5n[p{-,t)]/5p{x,t) and the 
fact that p{x,t) = for a; G dfl, which follows because 
the particles are confined to the box, has been utilized. 
We then restrict the perturbation field so as to satisfy 
4>{x,T ~ t) = (p{x,t). This allows us to compare ([15]) 
with (m , and we thereby obtain the following expression: 

TRx{x,T/2;x',t') ^ Cx,-v-j{x,T/2;x' ,t') 

+Cx^v-pG{x,T/2-x\t')/j, (16) 



where tilde denotes the time-reversal symmetric part, de- 
fined as, e.g., B{x,t;x' ,t') = {B{x,t;x',t') + B{x,t — 
t; x' , T — t'))/2. Next, setting e = 0, a similar calculation 
for (X{Xi t))o A under the assumption A{x, t) — A{x^ r — 
t) yields the following expression for Qy (a;, i; a;', t'): 

TQy{x,t/2-x'X) = tiY,-Vxj{x,T/2-x'X) 

+ tiY,vxpG{x,T/2;x',t')h- (17) 

Here, the equality Acr{x, t) = —Vxj{x, t), which follows 
from (fTTI) and (fT3|) . has been used. Equations (fT6|) and 
(fT7|l arc valid for arbitrary functionals, X and Y, and 
they express the relations among the response functions 
and correlation functions. 

We can obtain several important relations from (|16p 
and (Ull). First, the FRRs ® and ^ can be obtained 
by setting F{x) — and preparing the system with the 
canonical distribution, Pdv] = exp{—PH[(p]) with 
Z = J Dip exp{—f3'H[(f]). In this case, if X and Y sat- 
isfy time-reversal antisymmetry, i.e., = — 
and l^[[/o]] = — where p{x, t) = p{x, t — t), it can 
be confirmed that the second terms on the r.h.s. of (fT6|) 
and P?)) vanish, as a consequence of the detailed bal- 
ance condition Po{M]\[Po])Pc[Po] = Pa{[[p]]\[Pr])Pc[PTl 
where p^ = p{-,t). Because V[[p]](a;,t) and W[[p]](a;, t) 
satisfy time-reversal antisymmetry at i = t/2, we obtain 
dH) and © from ^ and ^ by setting X = V, Y = W 
and t — t/2. 

In contrast to the case analyzed above, when F{x) ^ 
and/or Pi ^ Pc, the system is not in equilibrium. In this 
case, we can obtain ([5]) from (fT6|) and (fT7|) as follows. 
First, if we set a;' = a; and t' = t/2, HI]) and ([17]) yield 

j[Cx.-v-j{x,t-x,t) -TRx{x,t;x,t)] = 

-{X{x,t)V-[p{x,t)G{x,t)])^, (18) 
jtr[DY.^~Vxj{x,t;x,t) -TQY(x,t;x,t)] = 

-{Yix,t)-[\7 X pix,t)G{x,t)]}„,{19) 

where t' has been replaced with t. We then integrate the 
sum of (fT51) and ITOl) over Q to obtain 



7 / dx{Cx.-v jix,t\ x,t) — TRx{x,t; x,t) 
Jn 

+tT[DY,^Vxj{x,t;x,t) - TQY{x,t;x,t)] } = 

l^j dxp{x,t)V{x,t)-G{x,t)^^, (20) 

where V = VX — VxY and the fact that p{x, i) = for 
X £ dfl has again been utilized. We next set X = V and 
y = VV. Because the relation W{x, i) - V x \V{x, t) = 
j{x,t)/ p{x,t) follows from ([6]) and ([7]), the integral on 
the r.h.s. of (l20l) can be rewritten as 



/ 

Jn 



dxj{x,t) ■ G{x,t) = Jn{t), 



(21) 
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where we have used the fact that j(x, t) = for x € dU. 
FinaUy, combining ((20)l and pT|) . we obtain the desired 
expression given in ([5]). 

In this way, we have obtained a concise relation be- 
tween the rate of energy dissipation and the magnitude 
of the FRR violation for field variables. Before closing 
this Rapid Communication, we make several remarks on 
the significance of this result. 

A model similar to ^ has also been studied in our 
previous report There, an equality relating the dissi- 
pation rate and the magnitude of the violation of the 
FRR for each mechanical degree of freedom was de- 
rived. The present result ([5|) contains the previous re- 
sult in Ref. [91] as a special case. If the spatial resolution 
of our measurement is sufficiently good so that we can 
distinguish every particle, the density field can be de- 
fined using the coordinates of the particles as p{x,t) = 
127=1 ^ Xi{t)). In this case, the entropy in ([72]) . i.e., 
— J^dxp{x,t)[\n p{x,t) — 1], becomes independent of the 
configuration of the particles due to the delta-function 
nature of the density field. Therefore, the contribution 
of the entropy term to the total dissipation (jH) vanishes, 
and the dissipation arises from the purely mechanical ori- 
gin. In this situation, we can prove the mathematical 
equivalence between the pres ent result ([5]) and the previ- 
ous result in (see 2^, 13]). In contrast, if our spatial 
resolution is restricted and we cannot distinguish each 
particle, the entropy changes in time and it measures 
the loss of information on the arrangement of the parti- 
cles. In this case, the variation of the entropy term also 
contributes to the dissipation ((4]). This is in accordance 
with Sekimoto's discussion that the definition of dissipa- 
tion should change depending on the level of description 
3l( . Our analysis here remains valid even in this case, 
and the FRR violation in the r.h.s. of ([5]) can capture the 
entropic component of the dissipation. 

The present result in ^ has also a practical advantage 
compared to the previous result in [9|]- In our previous 
result, in order to measure the FRR violation, one has 
to track each particle and to apply spatially uniform per- 
turbative force separately to every particle. This is quite 
difficult to be experimentally achieved. In contrast, mea- 
surement of the FRR violation in the present equality ^ 
does not require the identification of each particle. More- 
over, we can utilize spatially nonuniform fields, (f) and A, 
that can be applied to the entire system at one time. 
The latter should be more easily done in experiments. 
Therefore, the present approach extends the possibility 
to study the energetic aspect of nonequilibrium systems 
through measurement of the FRR violation. 
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